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ABSTRACT 

We report results from HST-NICMOS imaging of the extremely dense core of the globular 
cluster Terzan 5. This highly obscured bulge cluster contains a low-mass X-ray binary and at 
least two millisecond pulsars, with evidence for a large millisecond pulsar population. It has been 
estimated to have one of the highest collision rates of any galactic globular cluster, making its 
core a particularly conducive environment for the production of interacting binary systems. 

We have reconstructed high-resolution images of the central 19" x 19" region of Terzan 5 by 
application of the drizzle algorithm to dithered NIC2 images in the F110W, F187W, and F187N 
near-infrared filters. We have used a DAOPHOT/ALLSTAR analysis of these images to produce 
the deepest color-magnitude diagram (CMD) yet obtained for the core of Terzan 5. We have 
also analyzed the parallel 11" x 11" NIC1 field, centered 30" from the cluster center and imaged 
in F110W and F160W, and an additional NIC2 field that is immediately adjacent to the central 
field. 

This imaging results in a clean detection of the red-giant branch and horizontal branch in 
the central NIC2 field, and the detection of these plus the main-sequence turnoff and the upper 
main sequence in the NIC1 field. We have constructed an H versus J—H CMD for the NIC1 
field, which provides a measurement of the infrared reddening, E(J—H) = 0.72 (corresponding 
to E(B — V) = 2.16). We obtain a new distance estimate of 8.7 kpc, which places Terzan 5 within 
less than 1 kpc of the galactic center. The NIC1 CMD provides an apparent detection of the red 
giant branch bump, about 1 mag below the horizontal branch, indicating that the metallicity is 
at least solar and possibly somewhat higher. A number of blue objects are detected in the central 
NIC2 field and the NIC1 field, which hint at the presence of a blue horizontal branch. We have 
also determined a central surface-density profile which results in a maximum likelihood estimate 
of 7" 9 ± 0"6 for the cluster core radius. 

We discuss the implications of these results for the dynamical state of Terzan 5. 

Subject headings: globular clusters: individual (Terzan 5) — Hertzsprung- Russell diagram — infrared: stars 
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1. Introduction 

The extremely dense core of the bulge globular 
cluster Terzan 5 makes it a likely location for the 
dynamical production of interacting binaries. It 
has been estimated to have the highest stellar in- 
teraction rate of any galactic globular cluster (Ver- 
bunt & Hut 1987). Indeed, Terzan 5 is one of 
only 12 clusters known to contain a low-mass X- 
ray binary (LMXB; Johnston, Verbunt, & Hasinger 
1995); this object exhibits both transient behavior 
and bursts. Terzan 5 also contains an eclipsing mil- 
lisecond pulsar (MSP; Nice & Thorsett 1992), many 
steep-spectrum radio point sources within 30" of 
the cluster center, and steep-spectrum extended ra- 
dio emission in the cluster core (Fruchter & Goss 
1995). Lyne et al. (2000) have recently detected 
a second MSP just 10" from the cluster center. 
Based on the unresolved radio flux from the core 
of Terzan 5, Fruchter & Goss (2000) argue that it 
may have the largest pulsar population of any clus- 
ter in the Galaxy. They also detect three discrete 
radio sources that they interpret as likely pulsars. 

It has long been thought that dynamical interac- 
tions among stars in dense environments may lead 
to the production of interacting binaries and/or 
merged stars. Interaction mechanisms include two- 
body tidal capture and various three-body interac- 
tions, such as exchange reactions in which a mas- 
sive degenerate star may replace one member of a 
main-sequence binary. This picture is supported by 
Hubble Space Telescope (HST) studies, which have 
made the first detections of cataclysmic variables 
(CVs) in cluster cores (e.g. in NGC 6397; Cool et 
al. 1995) and have also detected substantial, cen- 
trally concentrated blue straggler populations (e.g. 
in M80; Ferraro et al. 1999a). 

Ground-based study of Terzan 5 has been ham- 
pered by its high obscuration (Ay ~ 7.7; Ortolani, 
Barbuy, & Bica 1996) and the extreme crowding 
in its core. Given its galactic coordinates, (£ = 
3°8,b = 1?7), the line of sight to Terzan 5 passes 
within 0.6 kpc of the galactic center. Ortolani et 
al. (1996) produced the first color-magnitude dia- 
gram (CMD) for Terzan 5, from imaging obtained 
under conditions of exceptional seeing (0"3 - 075) 
at the New Technology Telescope (NTT). Never- 
theless, their I vs. (V — I) CMD reaches only to 
just below the horizontal branch (HB) in the cen- 
tral 0.'5 radius about the cluster center. Moreover, 



their CMD is strongly affected by differential red- 
dening, with the red giant branch (RGB) spanning 
nearly a magnitude in color. Based on their CMD, 
Ortolani et al. (1996) estimate that the metallicity 
of Terzan 5 is approximately solar, making it one 
of the most metal rich of all globular clusters. This 
is consistent with the range of earlier estimates of 
[Fe/H] = +0.24, based on integrated infrared pho- 
tometry (Zinn 1985), [Fe/H] = -0.28, based on in- 
tegrated spectroscopy (Armandroff & Zinn 1988), 
and [Fe/H] = +0.07, also based on integrated spec- 
troscopy (Bica et al. 1998). Ortolani et al. (1996) 
also obtained a distance estimate of d = 5.6 kpc 
from their CMD, suggesting that Terzan 5 may be 
significantly closer than previous estimates, which 
generally lie in the range d = 7 — 9 kpc. Given 
its galactic coordinates, this range allows the pos- 
sibility that Terzan 5 is located very close to the 
galactic center. 

Terzan 5 is a natural target for HST-NICMOS 1 
imaging studies, both to better characterize the 
stellar content and structure of the central regions, 
and to search for additional evidence of the pres- 
ence of interacting binaries there. We report here 
results from NICMOS observations of the central 
region of Terzan 5. We describe the observations 
in §2, present the analysis in §3, and discuss our 
conclusions in §4. This paper concentrates on re- 
porting the analysis of the cluster CMD and central 
structure. A companion paper reports the results 
of a search for variable stars and the counterparts 
of the bright LMXB and the eclipsing MSP (Ed- 
monds et al. 2001). While the present paper was in 
the final stages of the refereeing process, an inde- 
pendent NICMOS study of the CMDs of several red- 
dened bulge clusters, including Terzan 5, appeared 
in print (Ortolani et al. 2001). Their study focuses 
on cluster age determination by isophotc fitting; we 
compare our results for the age of Terzan 5 with 
theirs in §4.4. 

2. Observations 

The primary targets for our NICMOS observa- 
tions were the central region of Terzan 5 and the 
MSP position 0'5 from the center (Edmonds et 
al. 2001). We imaged Terzan 5 with two primary 
pointings, one in which NIC2 was centered on the 

1 NICMOS = Near Infrared Camera and Multi Object Spec- 
trometer 
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cluster center and one in which NIC1 was centered 
on the MSP position. These observations produced 
parallel NIC1 observations at 30" from the cluster 
center and parallel NIC2 observations of a field ad- 
jacent to and slightly overlapping with the central 
field. Figure 1 shows the three fields used in the 
present study, viz. the central and adjacent NIC2 
fields, and the NIC1 field that was obtained in par- 
allel with the central field. The exposures used in 
this study include: (1) eight exposures of the cen- 
tral NIC2 field in each of the filters F110W, F187W, 
and F187N, (2) eight exposures of the NIC1 field 
in each of the filters F110W, F160W, and F187N, 
and (3) two exposures of the offset NIC2 field in 
F187W. For the exposures of the NIC1 and cen- 
tral NIC2 fields, the durations were 320 s for the 
broad-band (W) filters and 1536 s for the narrow- 
band (N) filter. For the offset NIC2 field, the total 
exposure time was 832 s. 
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Fig. 1 — Map of the three NICMOS fields used in this 
study. North is up and east is to the left. The star 
indicates our determined center of Terzan 5 and the 
circle indicates our maximum-likelihood determination 
of the core radius. The two eastern fields were observed 
with the NIC2 (field width = 19'.'2), while the smaller 
western field was observed with the NIC1 (field width 
= 11"). 

The eight observations of the central NIC2 field 
and the parallel NIC1 field were dithered, in or- 
der to increase the effective angular resolution, and 
consequently the photometric accuracy, in this ex- 



tremely crowded field. The F187W and F187N 
filters, which are centered on Paschen-a, were 
chosen for identifying emission-line objects. The 
F110W filter (which approximately corresponds to 
the Johnson J band) was included to produce a 
broad-band CMD, in combination with the F187W 
filter. Since F187W is not available for NIC1, we 
instead used F160W (which approximately corre- 
sponds to the Johnson H band) for the redder 
broad-band filter for this camera. While this has 
the disadvantage of providing a somewhat different 
CMD for the NIC1 field, a substantial advantage is 
that it allowed us to construct a NICl-field CMD 
in a standard ground-based filter system (J and 
H). This proved to be quite valuable in allow- 
ing us to determine new values for the metallicity, 
reddening, and distance of Terzan 5, as described 
below. 

The eight dither positions for the NIC2 field 
were selected in a 3 x 3 array (minus one cor- 
ner position). The basic step between dither po- 
sitions was 3 NIC2 pixels (0"23), with a random 
0.125 NIC2 pixel (0"009) step added to provide 
subpixel sampling. Each dither position was ob- 
served for one HST orbit. The 8-orbit observation 
sequence was split over two visits, separated in time 
by one month. Cross correlations of the resulting 
images indicate that subpixel pointing was accu- 
rately maintained from orbit to orbit within each 
visit. However, there were pointing shifts by about 
0.5 and 0.2 NIC2 pixels in the two orthogonal di- 
rections between the visits. These additional shifts 
acted in concert with the intentional 0.125 pixel 
shifts to provide subpixel sampling. 

3. Data Analysis 

3.1. Drizzle Reconstruction of Images 

The extraordinarily crowded central NIC2 field 
presents a substantial challenge to photometric 
analysis. The NIC1 and NIC2 provide full Nyquist 
sampling (i.e. 2 pixels per resolution element) down 
to l,um and 1.7/xm, respectively. Thus the NIC2 
images in the F110W band are strongly undcrsam- 
pled. The effect of this undersampling is alleviated 
by our dithering strategy for the central region. 
In order to make optimal use of the information 
contained in the dithered images, we employed the 
"drizzle" (variable-pixel linear reconstruction) al- 
gorithm of Hook & Fruchter (1997) to recover the 
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full resolution provided by the HST optics. We 
used the implementation of this algorithm from 
the STSDAS extensions of IRAF 2 , as described by 
Hook, Pirzkal, & Fruchter (1999). 

For the reconstructed NIC1 and NIC2 images, 
we chose pixel sizes of 0'.'022 and 0'.'038, respec- 
tively, which are half of the original pixel sizes. We 
chose a value of 0.5 input pix for the "drop" size; 
this parameter determines the magnification of the 
projection of each input pixel into the output im- 
age. Examination of the output drizzled images in- 
dicated a moderate level of residual pixel-to-pixel 
fluctuation, suggesting incomplete PSF sampling. 
These fluctuations were particularly evident in the 
vicinity of stellar image cores in the residual image 
on performing PSF-subtraction photometry. We 
reduced this noise by applying a Gaussian smooth- 
ing filter (<r = 0.75 pix) to the drizzled images. 
Figure 2 shows a comparison of a single F110W ex- 
posure of the central region and the result of com- 
bining the eight F110W images using the drizzle 
algorithm. 

Fig. 2. — Comparison of the central 5" x 5" section of a 
single F110W NIC2 exposure of the cluster center (left 
panel) with the corresponding section of the image that 
results from a drizzle reconstruction from eight dithered 
images (right panel). 

3.2. Photometry 

Given the substantial increase in the resolu- 
tion of the PSF that results from the application 
of the drizzle algorithm, it is feasible to carry 
out PSF-subtraction photometry We employed 



2 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



the DAOPHOT/ALLSTAR software package (Stet- 
son 1987; Stetson & Harris 1988) to determine an 
empirical PSF for each frame and simultaneously fit 
this to all detected stars. Approximately 100-150 
stars were used to determine a PSF with quadratic 
spatial variation across the frame. A complication 
in the process lies in distinguishing between faint 
stars and artifacts of the complex PSF structure of 
bright stars. We used visual inspection of apparent 
detections in the vicinity of bright stars to help fil- 
ter out artifacts. Since the diffraction pattern has 
a strong wavelength dependence, inter-comparison 
of the F110W and F160W images proved useful in 
this filtering process. A PSF radius of 20 pixels 
was used for both the NIC1 and NIC2 fields. This 
size was chosen as a compromise between the goals 
of subtracting as much of the halos of bright stars 
as possible, while still retaining a sufficient num- 
ber of stars, with no brighter neighbors within the 
PSF radius, to accurately characterize the spatial 
variation of the PSF. 

Given the higher degree of crowding in the cen- 
tral NIC2 field, relative to the offset NIC1 field, 
and the nearly twice as large NIC2 pixels, PSF 
subtraction was less effective in the NIC2 field. A 
clear indication of this is that the limiting magni- 
tude is considerably fainter in the NIC1 field for 
the same amount of exposure (see §4). Also, ex- 
amination of the residual frames for the NIC1 and 
NIC2 fields indicates a considerably greater amount 
of residual light in the latter. Nevertheless, PSF- 
subtraction photometry produced significantly bet- 
ter results than simple aperture photometry, for 
both the NIC1 and NIC2 frames. 

3.3. Calibration 

The NICMOS pipeline calibration system pro- 
duces images that are calibrated in flux units, i.e. 
the images are given in count rate (DN s -1 ). These 
can be converted to physical flux units (Janskys or 
erg s _1 cm~ 2 A -1 ) by multiplying by a correspond- 
ing scale factor that is tabulated for each filter. 
These fluxes can then be expressed in magnitude 
units in the AB„ or Vega systems; we use the lat- 
ter. As noted previously, the F110W and F160W 
filters, which were used for the NIC1 field, are 
analogs to the ground-based J and H bands. We 
computed a linear transformation between the HST 
and ground-based systems using NICMOS photom- 
etry for five standard stars with a wide color range 
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(-0.1 < J-H < +2.1) that are listed in the NIC- 
MOS Data Handbook, Version 4.0. Preliminary 
NICMOS photometry for the stars is given at the 
NICMOS website. The resulting transformations 
are: 

J = F110W + 0.033 - 0.335 (FU0W - FW0W) 

(1) 

H = FimW - 0.001 - 0.092 {F11QW - F160W) 

(2) 

These results are almost identical to those of 
Schulte-Ladbeck et al. (1999) who followed a sim- 
ilar calibration procedure. The RMS residual is 
0.16 mag for the J-transformation and 0.05 mag for 
the ^-transformation. The large value of the for- 
mer is mostly determined by one large residual (for 
the red standard BRI0021 with J - H = 1.17); the 
RMS of the other four J- residuals is 0.07 mag. We 
also note the existence of a significant color term 
in the J-transformation, which also contributes to 
the overall uncertainty. 

The transformations given by equations (1) and 
(2) allowed us to construct a standard near-infrared 
CMD for the NIC1 field. Since F187W (rather 
than F160W) serves as the redder broad-band fil- 
ter in the central NIC2 field, it is necessary to use a 
color index for this field that does not have a strict 
ground-based analog. 

4. The Color-Magnitude Diagram 

In constructing CMDs, we required that the cen- 
troids of stars determined by ALLSTAR match to 
0.5 pix between the F110W and F187W frames. 
This tight limit on the positional offset helped to 
filter out both weakly detected stars with uncertain 
photometry and PSF artifacts. We also eliminated 
stars that have neighbors with comparable mag- 
nitudes and centroids within 3 pixels, since these 
were strongly blended and thus produced substan- 
tial photometric uncertainty. 

Figure 3 shows a comparison of the CMDs ob- 
tained for the NIC2 and NIC1 fields. The ordinate 
is F110W in both cases. It is immediately obvi- 
ous that the central NIC2 field is much more popu- 
lous, owing to the much higher stellar density and 
the four times larger area. Not surprisingly, the 
offset NIC1 field produces much deeper photome- 
try, given the lower crowding and the smaller pixel 
size. The horizontal branch (HB) is clearly detected 
in both fields at a magnitude of F110W w 15.9, 



demonstrating the consistency of the calibration 
between the NIC1 and NIC2 cameras. The NIC1 
photometry also clearly detects the main sequence 
turnoff (MSTO) at F110W w 20 and upper main 
sequence; the faint-end limit of the NIC2 photom- 
etry is brighter than this at F110W w 19. 
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Fig. 3. — Comparison of the CMDs for the central 
NIC2 field (left panel) and the NIC1 field (right panel). 
The ordinate is F110W in both cases, while the abscissa 
is the available broad-band color index. The much 
greater depth of the NIC1 CMD is a consequence of the 
much lower crowding and 2x higher resolution in this 
field. The circled objects in the NIC2 field are photo- 
metrically variable (Edmonds et al. 2001); the brighter 
object is an RRab Lyrae and the fainter one is likely 
an eclipsing blue straggler. 

While one of the original goals of this study was 
to reach below the MSTO in the central NIC2 field, 
in order to survey the core of Terzan 5 for CVs, 
this is clearly not possible. Even with extensive 
dithering, the extraordinary crowding of this region 
substantially limits the accuracy of the photome- 
try, particularly for F110W > 18, which is just two 
magnitudes below the HB. Nevertheless, the pho- 
tometry in the central NIC2 field substantially su- 
persedes the quality of the excellent ground-based 
V and I photometry reported by Ortolani et al. 
(1996). Examination of their CMD for the cen- 
tral 26" radius region indicates that the limiting 
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magnitude is not much more than one magnitude 
below the HB and that the width of the RGB is 
A(V— I) = 0.92; they ascribe this to differential 
reddening. The greatly increased depth and pre- 
cision of the CMD that is afforded by our NIC- 
MOS observations permits us to search for unusual 
blue objects below the HB, if not below the MSTO. 
Furthermore, as has been recently demonstrated by 
Davidge (2000, 2001), near-infrared CMDs provide 
valuable information on the mctallicities, redden- 
ings, and distances of low-latitude globular clusters. 

In addition to constructing broad-band CMDs, 
we also constructed a F187N — F187W versus 
F187W diagram to search for evidence of emission- 
line objects. We find no convincing detections of 
such objects in either the central NIC2 field or the 
NIC1 field, in agreement with the analysis of the 
same data set by Edmonds et al. (2001, see their 
Fig. 8b). 

4.1. Blue Stars below the HB 

Figure 3 indicates a number of stars that lie blue- 
ward of the RGB in the NIC2 field. Edmonds et al. 
(2001) have shown that two of these blue stars show 
significant variability; these are circled in the left 
panel of Figure 3. As discussed by Edmonds et al. 
(2001), the brighter of these variables is an RRab 
Lyrae, while the fainter is likely an eclipsing blue 
straggler or possibly the LMXB counterpart. These 
two stars are both located at 8" from the cluster 
center. This puts the objects at the projected core 
radius. While the possibility of a chance superpo- 
sition always exists, particularly in a crowded field, 
the small radial offset of these stars from the clus- 
ter center makes a superposition of a nonmcmbcr 
unlikely. Moreover, given the consistency of the ap- 
parent magnitude of the RR Lyrae with the value 
expected for one located in Terzan 5, it appears 
likely that the star is at the cluster distance. While 
the photometry of the RR Lyrae appears fairly se- 
cure, the candidate blue straggler is clearly in the 
region where photometric uncertainty substantially 
widens the RGB. No variability is observed from 
the other two stars that lie close to the RR Lyrae 
in the CMD. 

Given the extreme crowding of the central field 
of Terzan 5, imperfect PSF-subtraction photometry 
can produce large photometric error as well as out- 
right artifacts, e.g. PSF "tendrils" that are detected 
as stars. We performed several tests of the photom- 



etry of the apparently blue stars, including visual 
inspection of the detections in the drizzled frames, 
variations in the photometric procedure, and cuts 
on the magnitude error estimates provided by ALL- 
STAR. Visual inspections generally proved useful at 
and above the HB level, where photometric prob- 
lems were fairly clear. However, by about 1 mag 
below the HB, it became quite difficult to discern 
photometric problems in this way, due to the ex- 
treme crowding of the field and the complexity of 
the PSF. Thus, we deemed automated photometric 
tests to be more useful for faint stars. To this end, 
we investigated cuts on the formal magnitude error 
estimates provided by ALLSTAR. The result of one 
such experiment is shown in Figure 4. In this case, 
a limit of 0.03 was placed on the formal errors of the 
F110W and F187W magnitudes. The three bluest 
stars are most likely foreground objects; a strong 
sequence of such objects is seen in the 2'.2 x 2'.2 
field surveyed by Ortolani et al. (1996). 
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Fig. 4.— CMD for the central NIC2 field. The pho- 
tometry is the same as in the left panel of Fig. 3, with a 
0.03 mag limit placed on the formal photometric error. 
The circled objects are the variables noted in Fig. 3; 
the fainter one does not pass the magnitude error filter. 
The dashed line indicates the location of the possible 
blue horizontal branch. The three bluest stars, near 
F110W = 17 are probably foreground stars. 



6 



Figure 4 provides a hint of the presence of a 
blue horizontal branch (BHB) that is parallel to 
the RGB and about 0.5 mag bluer in color. This is 
the expected location of the BHB in infrared col- 
ors, as can be seen in the photometry of four metal- 
poor clusters by Davidge & Courteau (1999). The 
hint of a BHB in Tcrzan 5 is intriguing, since the 
presence of a BHB in an extremely metal rich clus- 
ter is unusual. Rich et al. (1997) first discovered 
substantial BHB populations in old, metal-rich sys- 
tems. They found a blue extension of the strong red 
horizontal branch (RHB) in the globular clusters 
NGC 6388 and NGC 6441, based on HST WFPC2 
BV imaging. These clusters have high [Fe/H] val- 
ues of —0.60 and —0.53, respectively. Both clus- 
ters also have extremely high central densities, with 
logpo (A4© pc~ 3 ) values of 5.3 and 5.4, respectively 
(Djorgovski 1993). Rich ct al. (1997) suggest that 
the presence of a BHB in these clusters may be 
the result of enhanced mass loss from giants in- 
duced by stellar interactions. However, they note 
that the BHB/RHB ratio does not show the radial 
gradient that might be expected if BHB stars were 
produced by dynamical interactions in the cluster 
cores. Moreover, Bedin et al. (2000) have shown 
that the extended BHB population in the cluster 
NGC 2808 is no more centrally concentrated than 
the RGB population. A complicating factor is that 
once a star loses a substantial amount of envelope 
mass, by either stellar interactions or normal giant 
winds, it will tend to move to larger orbital radius 
as a consequence of mass segregation. Given HB 
lifetimes of ~ 10 s yr, this mass segregation pro- 
cess has enough time to at least partially reverse 
an initial central concentration of BHB stars. 

Since Terzan 5 has an even higher central 
density (logpo = 5.8; Djorgovski 1993) than do 
NGC 6388 and NGC 6441, dynamical interaction 
effects should be even more important for Terzan 5. 
Given the handful of potential BHB stars detected 
by our photometry, it is not feasible to examine 
the radial dependence of the BHB /RHB ratio. This 
will require wide-field infrared imaging of Terzan 5. 

4.2. The RGB Bump and Metallicity 

Stellar evolution theory predicts the presence of 
a RGB "bump," a pilcup that results from a pause 
in the ascent of the RGB that occurs when the hy- 
drogen burning shell extends past the point of last 
maximum depth of the convective envelope (see 



e.g. Iben 1968). Fusi Pecci et al. (1990) demon- 
strated that the location of the bump on the RGB 
is correlated with the cluster metallicity. Ferraro 
et al. (1999b) have measured the location of the 
RGB bump in 47 globular clusters and have cali- 
brated the relation between the offset of the bump 
magnitude from the HB magnitude (AV^™ P = 
Vbump — Vhb) and several different metallicity in- 
dices. The bump magnitude is brighter than the 
HB for mctallicities below [Fe/H] w — 1.2 and is 
fainter than the HB for higher metallicities. The 
clusters in the Ferraro et al. (1999b) survey range 
in [Fe/H] from -2.2 to -0.2. The high end of this 
range is of particular interest, for comparison with 
Terzan 5. The most metal rich cluster in the sam- 
ple is NGC 6528 with [Fe/H] = -0.23 on the Zinn 
(1985) scale and AV^ ip = +0.78. Davidge (2000) 
has recently reported near infrared photometry of 
NGC 6528, with the Canada-France-Hawaii Adap- 
tive Optics system, and finds a similar bump offset 
in the infrared, viz. 

Ferraro et al. (1999b) and Davidge (2000) note 
that examination of the cluster luminosity function 
provides a good means of determining the magni- 
tude of the RGB bump. The top and middle panels 
of Figure 5 show F110W luminosity functions for 
NIC2 and NIC1 fields, respectively. The results 
for the NIC2 frame are shown for both the entire 
frame and the corners of the frame only. The mo- 
tivation for separately considering the corners is to 
reduce the influence of photometric error, which 
is greatest in the central part of the frame owing 
to the greater crowding there. The bottom panel 
shows the ratio of the NIC2 and NIC1 luminosity 
functions, which provides a measure of the com- 
pleteness of the counts in the NIC2 field. Only 
bins containing at least three stars are shown. The 
solid horizontal line shows the overall count ratio 
for F110W < 18.4. The NIC2/NIC1 count ratio 
shows that significant incompleteness sets in in the 
NIC2 field for F110W > 18.0; this is about 2 mag 
fainter than the HB. This incompleteness is not 
surprising, given the extraordinarily high central 
density of Terzan 5, together with the complex NIC- 
MOS PSF structure. 
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Fig. 5. — Comparison of the luminosity functions for 
the central NIC2 field (top panel) and the NIC1 field 
(middle panel). The two histograms for the N1C2 field 
represent the results for the entire frame (upper his- 
togram) and the corners of the frame only (lower his- 
togram). Here the "corners" are defined as the region 
outside of an inscribed circle of radius 9'.'6. The bottom 
panel shows the ratio of the luminosity functions for the 
entire NIC2 and NIC1 fields. The solid horizontal line 
is the overall count ratio for F110W < f8.4 and the 
dotted horizontal lines represent the 1-a errors on the 
ratio. 

The HB is clearly evident at F110W = 15.9 in 
all of the histograms in Figure 5. There is also a 
smaller second peak at F110W = 16.9 evident in 
the NIC1 histogram. This feature corresponds to 
the clump of stars seen below the HB in Figure 3, 
and appears to represent the RGB bump. To es- 
timate its statistical significance, we compare its 
height (11 stars) to the mean of the two adjacent 
bins (3 and 4 stars); the result is a 2.1a t-ratio. 
Given that this is approximately where the RGB 
bump should occur in a very high metallicity clus- 
ter, we are fairly confident that we have detected 
it in the NIC1 field, fn contrast, examination of 
the histogram for the entire NIC2 frame shows no 
evidence of a RGB bump. There is a hint of a RGB 
bump in the histogram for stars in the frame cor- 
ners, although this feature has low statistical signif- 
icance (1.6cr). It may be that a small RGB bump in 



the NIC2 field has been smeared out by photomet- 
ric error. The lack of a RGB bump in the histogram 
for the entire NIC2 field is reflected by the dip in 
the NIC2/NIC1 count ratio at F110W = 16.9. 

The RGB bump detected in the NIC1 field pro- 
vides an estimate of the metallicity of Terzan 5. 
The middle panel of Figure 5 indicates a value of 
AFllOWy™ 15 = +1.0. We expect a similar value 
to hold for Ayg™ p , since there is little color differ- 
ence between the RGB bump and the intersection 
of the HB with the RGB. As noted above, the 
largest AV~Hg np observed by Ferraro et al. (1999b) 
is +0.78 for NGC 6528. For AV^ P = +1.0, their 
A^ u B mp - [Fc/H] Z85 calibration gives [Fe/H] Z85 - 
+0.26 for Terzan 5 (where Z85 = Zinn 1985). This 
is close agreement with the value of +0.24 from 
Zinn (1985), which is based on integrated infrared 
photometry. It is also consistent with the conclu- 
sion by Ortolani et al. (1996), based on the RGB 
curvature, that the metallicity of Terzan 5 is higher 
than that of NGC 6553 ([Fe/H] Z85 = -0.29) 3 and is 
probably solar. Thus there appears to be a strong 
indication that the metallicity of Terzan 5 is at least 
solar and is perhaps even somewhat higher. This is 
particularly interesting, in light of our detection of 
an RR Lyrae in the core (Edmonds et al. 2001, see 
§4.1 above), since this would make Terzan 5 the 
highest metallicity cluster in which an RR Lyrae 
has been observed. This appears to provide yet an- 
other indication of a high stellar interaction rate in 
the core of Terzan 5, on the assumption that colli- 
sionally induced mass loss was responsible for the 
production of the RR Lyrae. 

4.3. Reddening and Distance Modulus 

Figure 6 shows the CMD derived from the NIC1 
frame photometry, calibrated in the Johnson JH 
magnitude system. The upper RGB is sparsely 
populated, given the small size of this field and its 
30" offset from the cluster center. Nevertheless, the 
six stars above the HB are sufficient to character- 
ize the location and slope of the upper RGB. The 
dashed line represents a reddened fit of the Baade's 
Window (BW) M-giant sequence from Frogel & 



3 Several recent determinations of the metallicity of NGC 6553 
have been carried out by Barbuy et al. (1999), Cohen et al. 
(1999), and Coclho ct al. (2001). While there some differ- 
ences among these studies in detail, there is general agree- 
ment that the metallicity of NGC 6553 is near solar. 
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Whitford (1987, Table 3B), based on the Rieke & 
Lcbolsky (1985) infrared reddening law. The in- 
ferred infrared reddening is E(J—H) = 0.72 and 
the total infrared extinction is Ah = 1.18. This 
corresponds to a visual reddening of E(B — V) = 
2.18 and a F-band extinction of Ay = 6.75. Our 
reddening estimate is somewhat smaller than the 
value of E(B-V) = 2.49 obtained by Ortolani et 
al. (1996) from the position of the RGB in their 
ground-based VI CMD. Other previous redden- 
ing estimates include E{B — V) = 1.65 (Armandroff 
& Zinn 1988, based on a spectroscopic measure 
of an interstellar band), E(B-V) = 2.1 (Malkan 
1982, based on integrated infrared photometry), 
and E(B-V) = 2.42 (Djorgovski 1993, method un- 
specified). Inasmuch as the properties of Terzan 5 
are best determined in the infrared, it is interest- 
ing to note the close agreement between our new 
reddening estimate and that of Malkan (1982). 
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Fig. 6.— CMD for the offset NIC1 field. The field cen- 
ter is located at 30" from the cluster center. The pho- 
tometry has been transformed to the JH system. The 
dashed line is a fit of the Baade's Window M-giant se- 
quence, with E(J — H) = 0.72. The heavy solid curve is 
the CMD median ridge line and the lighter solid curves 
represent the first and third quartiles of the color dis- 
tribution. 

As we note in §1, Ortolani et al. (1996) obtained 
a significantly closer distance for Terzan 5 (5.6 kpc) 



than have other investigators. If we simply replace 
their reddening value of E{B — V) = 2.49 by our 
value of E{B — V) = 2.18, their distance estimate 
would be increased to 8.7 kpc. Alternatively, we 
can estimate the distance by comparison with the 
high-metallicity cluster NGC 6528, which has been 
studied in the infrared by Davidge (2000). The 
reddening-corrected HB magnitude of NGC 6528 
is Hub = 13.1 versus 13.0 for Terzan 5. Adopting 
a distance of 9.1 kpc for NGC 6528 (Harris 1996) 
also yields d = 8.7 kpc for Terzan 5. Recent galac- 
tic center distance determinations have produced 
values of dec = 8.3 ± 1.0 kpc (Carney et al. 1995) 
and d GC = 7.9 ± 0.3 kpc (McNamara et al. 2000). 
Since the line of sight to Terzan 5 passes within 
0.6 kpc of the galactic center, the total distance of 
Terzan 5 from the galactic center appears to be no 
more than about 1 kpc. Uncertainties in both the 
NICMOS calibration and the infrared reddening law 
contribute to the error in the distance determina- 
tion. These uncertainties are not well quantified. 
We note that if the distance modulus is accurate 
to at least 0.4 mag, then the uncertainty in the 
distance is less than about 2 kpc. 

The apparently small distance of Terzan 5 from 
the galactic center is interesting in light of the 
metallicity gradient in the bulge predicted by dissi- 
pative collapse models for its formation (e.g. Molla, 
Ferrini & Gozzi 2000). Frogcl, Ticde, & Kuchin- 
ski (1999) have measured a metallicity gradient of 
about —0.45 dex/kpe along the minor axis of the 
bulge. Feltzing & Gilmore (2000) have recently 
found a steeper gradient value of —1.3 dex/kpe. 
Based on near-infrared spectroscopy of bulge gi- 
ants, Ramirez et al. (2000) argue that there is lit- 
tle evidence for a gradient inside of 560 pc from 
the galactic center and find a mean metallicity of 
—0.21 ± 0.30 (s.d.) for the inner bulge. Since the 
metallicity that we find here for Terzan 5 is at least 
this large and the bulge metallicity drops outside of 
the inner bulge, an inner bulge origin for Terzan 5 
appears to be a reasonable inference. 

4.4. The Main Sequence Turnoff and Clus- 
ter Age 

As we noted in §4, it is evident from Fig- 
ure 3 (right panel) that the NICl-field photometry 
reaches below the main sequence turnoff (MSTO). 
The narrow 1.0 mag color range used in the NIC1 
JH CMD shown in Figure 6 emphasizes the spread 
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of the points about the mean color-magnitude re- 
lation. To determine the location of the MSTO, 
we computed the median color for 0.5 mag bins in 
H. We also computed the first and third quartiles 
of the color distributions in each magnitude bin, 
in order to characterize the width of the distribu- 
tions. These three color-distribution statistics are 
plotted as continuous curves in Figure 6. The me- 
dian color reaches its bluest value near H — 18.5, 
which is 4.3 mag below the HB. We identify this as 
the MSTO and note that its detection at 30" from 
the cluster center is an impressive demonstration of 
the capability of NIC1 imaging. The correspond- 
ing unreddened magnitude and color of the turnoff 
are H TO = 17.3 and (J-H) TO = 0.4. The semi- 
interquartile range of the color distribution at the 
turnoff is just 0.06 mag, indicating a fairly narrow 
photometric spread for a typical star. 

The magnitude difference between the HB and 
the MSTO, AV^g, is an important indicator of 
cluster age, as discussed by Chaboyer, Demarque, 
& Sarajedini (1996). Girardi et al. (2000) have 
recently computed evolutionary tracks for low- to 
intermediate-mass stars with a wide metallicity 
range, and have derived isochrones in the Johnson- 
Cousins system, including the J and H bands. 
These can, in principle, be used to calibrate the 
relation between AiJgg and cluster age. In ap- 
plications of the Ay^jg method of age determina- 
tion, it is usual to adopt a fixed absolute magni- 
tude for the HB. This is typically the RR Lyrae 
magnitude, which is assumed to be independent of 
age over the relevant age range for globular clusters 
(see e.g. Chaboyer, Demarque, & Sarajedini 1996; 
Rosenberg et al. 1999). For very high- metallicity 
clusters, where RR Lyraes may be deficient or ab- 
sent, an alternative choice for defining a fixed HB 
magnitude, is the zero-age HB (ZAHB) magnitude 
for a star of a chosen fudicial mass and appropriate 
metallicity (see Ortolani et al. 2001, as discussed 
below) . 

As a cluster evolves, the magnitude of the ZAHB 
is expected to become somewhat fainter, as indi- 
cated by the Girardi et al. (2000) isochrones, which 
are based on evolutionary tracks that include HB 
and post-HB evolution. The resulting evolution of 
the red giant clump (i.e. red HB) absolute mag- 
nitude is shown by Girardi & Salaris (2001, see 
their Fig. 1). Thus, an alternative procedure for 
calibrating the AiJgg-age relation is to use the 



evolving HB magnitude. In Figure 7, we compare 
these two calibration procedures for two values of 
the cluster metallicity. In the upper panel, we use 
the contemporaneous HB magnitude from the Gi- 
rardi et al. (2000) isochrones, while in the lower 
panel we use the fixed, zero-age HB magnitude for 
a star with an initial mass of 1M.Q, following Or- 
tolani et al. (2001). We define the contemporane- 
ous HB magnitude by the bluest point on the HB 
isochrones, since examination of the Girardi et al. 
(2000) isochrones indicates that stars linger near 
this location. It can be seen that the calibration 
curves flatten considerably for t > 10 Gyr when 
the evolving HB magnitude is used. In both pan- 
els, we show the uncertainty range produced by an 
estimated 0.2 mag precision in the determination 
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Fig. 7. — Magnitude difference between the horizon- 
tal branch and turnoff, AH^%, versus cluster age, for 
(a) a time-evolving HB magnitude, and (b) a fixed HB 
magnitude. In each panel, the curves are labeled by the 
[Fe/H] value. The solid horizontal line in each panel is 
the observed value of AJ/hb f° r Terzan 5 and the dot- 
ted lines represent the estimated range of uncertainty. 
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The curves in the upper and lower panels of Fig- 
ure 7 give similar results for the best-fit cluster age 
for a given value of the metallicity. However, it 
can be seen that the uncertainty ranges are signif- 
icantly different for the two calibration methods. 
For an evolving HB magnitude, the best-fit ages 
are 15 Gyr for solar metallicity and 10 Gyr for su- 
persolar metallicity, with lower limits of 8 Gyr and 
7 Gyr, respectively. The flattening of the calibra- 
tion curves in this case precludes placing upper lim- 
its on the ages for the estimated 0.2 mag precision 
of Af/^g. The limits are tightened significantly 
for a fixed HB magnitude, where the best-fit ages 
are 14+ 4 Gyr for solar metallicity and ll+ 3 Gyr for 
supersolar metallicity. 

We note that Ortolani et al. (2001) have inde- 
pendently determined an age for Terzan 5 by fitting 
isochrones to a CMD determined from their NIC1 
observations of a field 30" from the cluster center. 
They find an identical offset of 4.3 ± 0.2 mag be- 
tween the HB and the MSTO, at F110W, to what 
we find at H, from which they infer an age of 
14 ± 3 Gyr, using newly computed isochrones in 
F110W-F160W HST filter system. They assume a 
fixed ZAHB magnitude, taken to be that of a star of 
initial mass of 1 Mq- They also adopt solar metal- 
licity. Thus, they have treated the case represented 
by the lower curve in the lower panel of Fig. 7, for 
which our best-fit age is also 14 Gyr. However, their 
computed values of AFIIOWhb have a steeper de- 
pendence on age than do the values for A/i^B that 
we derived from the Girardi et al. (2000) isophotes, 
allowing Ortolani et al. (2001) to place somewhat 
tighter limits on the determined age (±3 Gyr versus 
±4.5 Gyr). This evidently results from a difference 
in the behavior of the isochrone sets. The good 
agreement between the best-fit values for the age 
of Terzan 5 (assuming solar metallicity), obtained 
here and by Ortolani et al. (2001), is encouraging. 
However, the assumption of a fixed HB magnitude, 
when not using RR Lyraes to define the HB, should 
be investigated further. 

5. Cluster Structure 

Trager, King, & Djorgovski (1995) have studied 
the spatial structure of Terzan 5 by determining 
the surface-brightness profile from an /-band CCD 
image. They determined a core radius value of 
r c = 10'.' 7 and a half-mass radius of r h = 55". Since 



small numbers of individual bright giants domi- 
nate the surface-brightness distribution, especially 
in the /-band, star-count profiles generally provide 
a better means of determining the cluster structure 
(see e.g. Sosin & King 1997). 

We have constructed a surface-density profile for 
Terzan 5, using star counts from the central and 
offset NIC 2 fields (see Figure 1). It is necessary to 
include the offset field, since the cluster core fills the 
central field. We used the F187W star counts, since 
this is the only broad-band filter common between 
these two fields. Since our exposures for the offset 
field arc not dithered, it was not feasible to per- 
form PSF-subtraction photometry for these frames. 
Instead, we used simple aperture photometry and 
computed a magnitude offset between the central 
and offset fields from the stars in the overlap re- 
gion. Since the depth of the star counts varies with 
distance from the cluster center, we adopted a uni- 
form limiting magnitude for stars to be included in 
the profile determination that corresponds to 2 mag 
below the mean HB, i.e. F187W = 15.9. Figure 5 
indicates that the central field is fairly complete to 
this limit. 

We used iterative centroiding, in a circular aper- 
ture with an initial radius equal to the half width 
of the NIC2 field (9'.'6), to determine the location of 
the cluster center. As the centroid shifted from the 
starting position at the frame center, we reduced 
the size of the aperture to keep it entirely within 
the central NIC2 field at each iteration. The iter- 
atively determined centroid was offset by Of! 7 from 
the frame center. Trials with a range of values of 
the limiting magnitude and aperture size indicate 
that the uncertainty in the location of the center 
is of order 1". Our determined center position 
corresponds to the equatorial coordinates (J2000) 
a = 17 h 48 m 04?8 and 5 = -24° 46' 45. Table 1 
compares this position with previous determina- 
tions. There is particularly good agreement be- 
tween the present center determination and those 
of Djorgovski & Meylan (1993, after correction; see 
Table 1 notes) and Fruchter & Goss (2000). The 
latter is the centroid of the unresolved radio flux 
from the core, which they suggest arises from an 
unresolved MSP population. The close agreement 
between this radio center and our present optical 
center is striking. 

We binned the star counts into a central disk 
of radius 0'.'75 and nine concentric, logarithmically 
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Table 1 

Comparison of Cluster Center Determinations. 



Study a (2000) 6 (2000) A (") 



Terzan (1971) 


17 


48 


04.6 


-24 


46 


50.9 


6.7 


Djorgovski & Meylan (1993) b 


17 


48 


04.9 


-24 


46 


45 


1.1 


Picard & Johnston (1995) 


17 


48 


04.6 


-24 


46 


48.5 


4.7 


Fruchter & Goss (2000) 


17 


48 


04.9 


-24 


46 


44.7 


1.3 


Present study 


17 


48 


04.8 


-24 


46 


45.0 





a Total offset relative to present study. 

b A misprint in the 5 minutes value has been corrected per private 
communication with S. Djorgovski. 

Note. — All positions are given with the same precision as in the 
original source. 



spaced circular annuli that extend to 29" from the 
cluster center; the annuli are incomplete beyond 9". 
The surface density in each annulus is computed by 
dividing the number of stars in that annulus by the 
area of the annulus that lies within the combined 
NIC 2 central and offset fields. Figure 8 shows the 
resulting surface-density profile. We have included 
the Trager et al. (1995) surface-brightness profile 
for comparison 4 , with a shift applied to it so that 
the two profiles agree near r = 10". It can be seen 
that both profiles shown in Figure 8 have relatively 
fiat cores, with a somewhat smaller core radius for 
our surface-density profile. 

In order to determine the core radius in a man- 
ner independent from the adopted binning, we em- 
ployed the maximum-likelihood method outlined 
by Sosin & King (1997). In this approach, a sur- 
face density profile function is fit to the star posi- 
tion data directly, without resort to radial binning. 
The likelihood for the observation of a particular 

4 The Trager et al. (1995) data set includes three spectral 
bands, B, R, and /. The I and R data cover the entire 
profile, while the B data only cover the outer halo. Since 
the I and R data each tightly define a profile that is slightly 
different from that defined by the other band, we chose to 
include only their B and / data in this comparison. 



set of stellar radial position values (r^) is given by: 

N 

L(a) = IJ/(ri;a) (3) 

»=i 

where / represents a parameterized form for the 
surface density profile (in units of stars per unit 
area) and a represents a set of parameters on which 
/ depends. It is convenient to work with logi, 
which converts the product in equation (3) to a 
sum. Rather than compute logL over a grid in 
parameter space, to locate the point of maximum 
likelihood, we found it more convenient to use an 
optimization routine. We used an implementation 
of Powell's method from Numerical Recipes (Press 
et al. 1992, §10.5), which does not require the com- 
putation of derivatives of the likelihood function 
with respect to the parameters. For estimating the 
confidence ranges for the parameter estimates, we 
used the "bootstrap" method, i.e. random resam- 
pling of the star list with replacement (Efron 1982; 
Press et al. 1992, §15.6). We implemented this by 
drawing 1000 random resamples and determining 
the maximum-likelihood values of the parameters 
for each resample. We verified that the distribu- 
tions of the parameter values determined from the 
set of resamples are closely approximated by Gaus- 
sians, with means equal to the best-fit parameter 
values for the original sample. The standard devia- 
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tions of these distributions provide robust estimates 
for the uncertainty of the best-fit parameter values. 
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Fig. 8. — Surface-density profile from F187W star 
counts (solid dots and error bars) compared with sur- 
face brightness profile from Trager et al. (1995, open 
diamonds). The solid line is a maximum likelihood fit 
of eqn. (4) to the unbinned F187W star counts. 

Following Lugger, Cohn, & Grindlay (1995), we 
used the following modihed-power-law form for the 
azimuthally symmetric surface density distribution, 



fir) = f 



1+ [ - 



a/2 



(4) 



Here r is a radial scale related to the core ra- 
dius and a is the power-law index of the outer 
part of the profile. The normalization parameter 
/o is determined by the condition that the inte- 
gral of f(r) over the area in which the stars are 
located remains constant as the values of the other 
two parameters are varied. For the usual definition 
of the core radius as the point at which the sur- 
face density decreases to half of its central value, 
r c = (2~ 2 / a — l) 1 / 2 r . A power-law slope of a = -2 
corresponds to an "analytic King model" (King 
1962). Lugger et al. (1995) have shown that equa- 
tion (4) with a = —1.8 provides an excellent fit to 
an actual King (1966) model with a central concen- 
tration of c = 2.0, for all but the outer halo. As a 



consistency check, we performed a least-squares fit 
of equation (4) to the Trager et al. (1995) profile. 
We find r c = 10'.'2 and a = -1.9. This core ra- 
dius value agrees with their value of 10'.'7 to within 
5%; the power-law slope is consistent with a King 
model. 

Our maximum likelihood fit of equation (4) to 
the NIC2 star counts produces r c = 7'.'9 ± 0'.'6 and 
a = —2.1 ±0.2. It can be seen in Figure 8 that this 
form provides an excellent fit. The core radius is 
about 25% smaller than that found by Trager et al. 
(1995), while the power-law slope is similar to that 
produced by our fit to their profile. Interestingly, 
in a preliminary report of their profile fits, Trager, 
Djorgovski, & King (1993) found a core radius of 
7'/9, in agreement with our value. The slope value 
of a « — 2 again indicates consistency with a King 
model. 
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Fig. 9. — Combined surface-density profile (as de- 
scribed in Fig. 8) with King model overlays. The curves 
correspond to the indicated c values from left to right. 
The core radius and central density are taken from the 
maximum likelihood fit shown in Fig. 8. Note that the 
best match is obtained for concentration parameter in 
the range 2.0 < c < 2.25. 

To further explore the limits on the central con- 
centration of Tcrzan 5 that can be obtained from 
model fitting to the combined profile shown in Fig- 
ure 8, we have overplotted a set of King models 
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with a range of c from 1.75 to 2.25, as shown in 
Figure 9. Rather than formally fitting the profiles, 
we simply adopted the core radius and central sur- 
face density values from our maximum-likelihood 
fits of the modified power-law profile. It can be seen 
in Figure 9 that the c = 2.0 King model provides 
the best overall fit to the other part of the profile as 
determined by Trager et al. (1995). Given the chal- 
lenge of determining the surface brightness profile 
at large distance from the cluster center, with the 
high background of the galactic bulge and variable 
extinction, we do not view this as a rigorous deter- 
mination of the central concentration of Terzan 5. 
Rather, we see it as a demonstration that a rela- 
tively high central concentration of c w 2 is cer- 
tainly plausible. 

For our adopted normalization of the profiles 
near r = 10", the central surface density is 60% 
higher for the NIC2 profile than for that of Trager 
et al. (1995). The combination of the higher sur- 
face density and the smaller core radius implies 
that the central density is about twice as large 
as the value of po = 10 5 76 .M pc~ 3 determined 
by Djorgovski (1993), i.e. p w 10 6 X Q pc~ 3 . 
This puts Terzan 5 at the top of the central den- 
sity distribution for all clusters, although the tab- 
ulated values for collapsed-core clusters are more 
properly viewed as lower limits. Following Djor- 
govski (1993), we may estimate the central relax- 
ation time of Terzan 5, by applying the scaling re- 
lation, i rc cx po°' 5 r^, t° his estimate. The result 
is t IC = 4 x 10 7 yr, which is at about the 30th 
percentile of the distribution for 124 clusters from 
Djorgovski (1993). Of the clusters with shorter cen- 
tral relaxation times than Terzan 5, only about one 
third have apparently "normal" cores, i.e. have a 
central structure that is well fit by a King model, 
like that of Terzan 5. Thus, Terzan 5 has both an 
extremely high central density and an unusually 
short relaxation time for a cluster with a noncol- 
lapsed core. 

6. Discussion 

Using NIC2 imaging of the central 19" x 19" re- 
gion of Terzan 5 and NIC1 imaging of a field 30" 
from the cluster center, we have constructed color- 
magnitude diagrams and produced new measure- 
ments of the cluster reddening, distance, and metal- 
licity. We have also determined the surface-density 



profile and performed a maximum-likelihood fit to 
measure the core radius and halo slope. Our results 
underscore the extraordinary nature of this cluster. 

The result that the reddening is somewhat 
smaller than found by Ortolani et al. (1996) places 
the cluster at a greater distance than the strikingly 
small value of 5.6 kpc that they found. This new 
reddening, together with our measurement of the 
J-band magnitude of the horizontal branch, place 
Terzan 5 at the 60% larger distance of 8.7 kpc and 
thus within less than 1 kpc of the galactic center. 
Our distance measurement is consistent with the 
range of 7 to 9 kpc produced by most studies prior 
to Ortolani et al. (1996). 

The probable detection of the red giant branch 
bump, about 1 mag below the horizontal branch, 
indicates that the metallicity of Terzan 5 is at least 
solar and possibly somewhat higher ([Fe/H]zs5 = 
+0.26). This places Terzan 5 at the very top of 
the metallicity distribution of all globular clusters. 
This suggests that Terzan 5 was formed in the 
vicinity of its present location near the galactic cen- 
ter. 

The detection of the main sequence turnoff has 
allowed us to estimate the cluster age using the 
magnitude difference between the turnoff and the 
horizontal branch. Our result for the best-fit age — 
14 Gyr for an assumed solar metallicity and a fixed 
HB magnitude — is the same as that recently ob- 
tained independently by Ortolani et al. (2001) from 
a separate NICMOS data set. However, our derived 
uncertainty in the age is somewhat larger than 
they obtained, owing to differences in the isochrone 
sets used. Even with their tighter age precision of 
±3 Gyr, it is only marginally possible to distinguish 
between "young" (~ 10 Gyr) and "old" (<~ 15 Gyr) 
ages for Terzan 5. We also note the role of metallic- 
ity in the age determination; Fig. 7 indicates that 
the inferred age may be reduced by 3-5 Gyr if su- 
persolar metallicity were adopted. Given the ex- 
traordinarily high metallicity of Terzan 5, a tighter 
limit on the cluster age would have important im- 
plications for the chemical history of the galactic 
center region, particularly if the age turned out to 
be fairly old. 

Improvements in the age estimate could be 
achieved by a combination of: (1) NIC1 mosaic 
imaging of the central region to improve the photo- 
metric precision and depth, (2) additional NICMOS 
imaging at about 0.5-1' from the cluster center, to 
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improve the statistics of the deep CMD, (3) NIC- 
MOS observations of more near-infrared standard 
stars, to improve the determination of the magni- 
tude system transformations, (4) additional study 
of infrared isochrones to better fit the observations 
of high-metallicity systems, and (5) a tighter limit 
on the metallicity of Terzan 5. 

Our maximum-likelihood fit of a modificd- 
power-law form to the surface-density profile within 
0.'5 of the cluster center has produced a new core 
radius measurement of r c = 7'.' 9, which is somewhat 
smaller than the previous determination by Trager 
ct al. (1995). Nevertheless, our results underscore 
the fact that Terzan 5 has a resolved core and is well 
fit by a King-model profile. Given the relatively 
short central relaxation time of t IC = 4 x 10 7 yr, it 
is interesting to ask whether Terzan 5 should have 
already undergone core collapse. For a single-mass 
cluster, the core collapse time is formally about 
300 times longer than the central relaxation time 
(Cohn 1980). This would suggest that Terzan 5 
may be just now at the brink of core collapse (see 
also Cohn & Hut 1984). Furthermore, Meylan & 
Heggie (1997) have argued that any cluster with 
central concentration c > 2 should be considered 
to be either on the brink of core collapse or just 
past it. 

While the concentration parameter of Terzan 5 
has not been directly determined, given the diffi- 
culty of determining the tidal radius of this highly 
obscured cluster, our comparisons with King mod- 
els suggest that c « 2. Given the estimated half- 
mass radius of Terzan 5 from Trager et al. (1995), 
the ratio of core to half-mass radius is r c /r^ = 0.15. 
Grabhorn et al. (1992) have shown that this ra- 
tio has a maximum value of about 0.01 - 0.04 in 
post-collapse dynamical models for clusters. This 
suggests that the core of Terzan 5 may still be 
supported in a pre-collapse state by energy re- 
lease ( "heating" ) from primordial binaries that in- 
teract dynamically with other stars. Alternatively, 
Terzan 5 may be in a post-collapse state with a 
"super-expanded" core, produced by the action of 
a residual population of primordial binaries. Such 
binaries were not considered by Grabhorn et al. 
(1992). 

If Terzan 5 is indeed near the point of core col- 
lapse, we might expect to see enhanced binary ac- 
tivity in its core. Ferraro et al. (1999a) have in- 
terpreted the huge blue straggler population that 



they detected in the high-density, noncollapsed core 
of M80 as evidence of vigorous binary heating, 
which is temporarily staving off imminent collapse. 
Specifically, they propose that the blue stragglers 
are the result of the hardening of both primor- 
dial binaries and binaries produced by tidal cap- 
ture and three-body interactions. The vast major- 
ity of blue stragglers present in Terzan 5 should 
lie in the central NIC2 field of our data set, given 
the extreme central concentration of blue strag- 
gler populations in dynamically evolved clusters. 
The likely eclipsing blue straggler identified by Ed- 
monds ct al. (2001) lies at about 8" from the clus- 
ter center. Highly dithered NIC1 photometry of 
the central region may aid in detecting additional 
blue stragglers in the core, although distinguish- 
ing between blue stragglers and photometric ar- 
tifacts in this crowded central region could prove 
challenging. Ground-based adaptive optics obser- 
vations with the Hokupa'a system on Gemini may 
well prove superior to NICMOS for probing the core 
of Terzan 5, given the large telescope aperture, the 
ft' 02 pixel size of the imager, the 20" field of view, 
and the ft/07 or better expected FWHM of the im- 
ages under good seeing conditions. 

Multi- wavelength studies of Terzan 5, including 
additional searches for millisecond pulsars and X- 
ray sources, will certainly be useful for assessing 
the level of stellar collisional activity in the core. 
Our group has recently obtained a deep Chandra 
observation of Terzan 5. However, the LMXB was 
in a high state during this observation, which may 
preclude the detection of a large population of faint 
X-ray sources in the central region. In any case, the 
X-ray spectral observations should provide a mea- 
surement of the hydrogen column density towards 
Terzan 5, and thus an independent estimate of the 
obscuration. This will provide a check on our new 
distance estimate. Our NICMOS study underscores 
the extraordinary physical conditions present in the 
core of Terzan 5. Thus, it remains of great interest 
to test the prediction that this region is a factory 
for the production of substantial numbers of inter- 
acting binaries. 
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